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The nonadiabatic electron-phonon corrections for the superconducting pairing are investigated 
for a specific tight-binding model corresponding to a 2d square lattice. This permits to investigate 
the role of various specific properties like the band filling, nesting effects and a realistic van Hove 
singularity on the superconducting effective pairing beyond Migdal's limit. The main results are: 
(i) Starting from a momentum independent electron-phonon coupling the nonadiabatic effects lead 
to an effective pairing which is strongly dependent on frequency and momentum, (ii) If instead 
the electron-phonon coupling is mainly forward (as due to correlation effects) the resulting pairing 
results to be strongly enhanced. These results confirm but also extend the simplified scheme used 
up to now to compute these properties. In this respect our results link the nonadiabatic effects to 
the specific properties of realistic materials. 
PACS rmmbers;74.20.Mn, 71 .38.-Hi, 63.20.Kr 



I. INTRODUCTION 



Unconventional superconductors, such as, cuprates and fuUerides, are characterized by Fermi energies Ep much 
smaller than those of conventional metals.ETcl In this situation, the energy scale loq of the mediator of the super- 
conducting pairing can be comparable to Ep so that the quantity ujq/Ef is no longer negligible as it happens for 
conventional superconductors. As a consequence, vertex correcjbions in the normal state electronic self-energy become 
relevant because of th(i.non- validity of Migdal's approximationH and also the Migdal-Eliashberg (ME) approach to the 
superconducting statecld is drastically modified.l3 

In recent years we have generaJizfid both the superconducting and the normal states in order to include the first 
corrections beyond Migdal's limit .Q^llil We considered the energy scale to be given by phonons so that the parameter 
luq/Ep measures the degree of adiabaticity of the lattice dynamics compared to the electronic kinetic energy. For 
luq/Ep ^ 1 and electron-phonon {e-ph) coupling A < 1 one recovers the ME theory of the electron-phonon coupled 
system. By enhancing the coupling A the electron develops a large effective mass (polaron) and eventually forms 
bi-polaronic states. However, the e-ph coupled system may display a regime which is different from both the ME 
picture and the polaronic (bi-polaronic) regime. We have named this situation nonadiabatic, and it is characterized 
by quasi- free electrons, i.e., A < 1, coupled in a nonadiabatic way {ojo/Ep not negligible) to the lattice. We have 
shown that in this regime it is possible to define a perturbative theory where the expansion parameter is roughly 
given by Xloq/Ep. Moreover, the corrections arising from the nonadiabatic hypothesis have a non-trivial structure in 
both the momentum transfer q and the Matsubara exchanged frequency uj. In fact, the vertex correction appearing 
in the normal state self-energy becomes positive (negative) when vpq < oj {vpq > oj), where vp is the Fermi velocity.u 
The generalization to the superconducting transition reveals that this situation is also encountered for the class of 
diagrams beyond Migdal's limit relevant for the Cooper channel. Concerning the critical temperature Tc, as long as 
the momentum transfer is less than loq /Vfri the nonadiabatic corrections lead to a strong enhancement of Tc also for 
moderate values of the e-ph coupling A.QlI The opening of a superconducting gap changes drastically the behavior 
in momentum-frequency space of the vertex especially in the low temperature limit. In this regime the vertex can be 
positive also when vpq > cu probably leading to anomalous values of the gap-to- Tc ratio.til 

The results presented in the previous papers on the nonadiabatic e-ph interaction theory have been derived by 
employing higly simplified approximations. Although several results were approximation- independent, it is important 
to study the behavior of the nonadiabatic corrections by using realistic models. The aim of this paper is to provide 
a study of both the vertex and cross diagrams in the framework of a specific tight-binding model for the electronic 
band. The advantage of such a study is twofold: it permits to confirm the results obtained by using simple models but 
also the whole picture is made much richer so that it is possible to relate specific nonadiabatic properties to specific 
features of the material. Finally, this paper can be also regarded as a preliminary work for a completely numerical 
solution of the Eliashberg equations generalized beyond Migdal's limit. 

The paper is organized as follows. In Sec.|| we introduce the equations for the superconducting instability driven 
by the e-ph coupling generalized in order to include the first nonadiabatic corrections represented by the vertex and 
cross diagrams for a specific tight-binding model. In Sec. Ill wc study the vertex and cross corrections for a model 
in which the e-ph interaction is momentum independent while in Sec.|^ we analyze the effect on the nonadiabatic 
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corrections of an e-ph coupling peaked at small values of the momentum transfer as induced by electronic correlations. 
The last section is devoted to the conclusions. 



II. BEYOND MIGDAL-ELIASHBERG THEORY 



The presence of nonadiabatic corrections beyond Migdal's limit affects both the normal and superconducting prop- 
erties of the metallic systems in a non-trivial way. In fact, as it has been already shown in Refs. up to the first 
corrections of order Xloq/Ep, the normal state electronic self-energy is modified by only one nonadiabatic correction, 
the vertex correction, while for the superconductive instability one needs to include also the cross correction. This 
situation is shown in Fig.|l| where the normal state electronic self-energy S^r (Fig.|l|a) and the linearized gap equation 
(Fig.|l|b) are depicted in order to include the first corrections beyond Migdal's limit. 

In order to find the critical temperature Tc, one has to solve the following coupled set of equations relating the 
diagonal, I]Ar(fc), and the off-diagonal, Ss(fc), electron self-energies: 

SAr(fc) = ^o{k - k')[l + P{k, k')]G{k'), (1) 
k' 

Ss(fc) = J2 fc')G(fc')G(-fc')Ss(fc'), (2) 

L./ 



where k = {iivk, k), k' = {iujk' , k'), and we use the notation J2k — J^l^^ Sk- above equations, Vo{k — k') = 

g(k — 'k')'^D{k — k') where (^(k — k') is the e-ph matrix element, G{k) and D{q) are the electron and phonon propagators 
given by: 

G-\k)=iuJk~e{k)-j:N{k), (3) 



and 



D{q) = (4) 



-I 



respectively. In equations (j^) and (|4|), e(k) is the electronic band dispersion and wq is the phonon frequency which 
we assume dispersionless for simplicity. In Eq.(|l|), P{k, k') is the first vertex correction beyond Migdal's limit, which 
according to Fig.|l|a is given by: 

P{k, k') = ^oik - p)G{p -k + k')G{p). (5) 
p 

The total pairing interaction V{k^ k') appearing in the linearized gap equation is instead given by: 

V{k, k') = V^ik - k')[\ + 2P{k, k')] + C(/c, k'), (6) 
where P(fc, k') is again the vertex correction d) while C(/c, k') is the cross term (last diagram appearing in Fig j^b) : 

C{k, k') = YMk- p)Vo{p - k')G{p)G{p ~k- k'). (7) 
p 

The ordinary Migdal-Eliashberg equations are obtained by neglecting both the vertex P and cross C corrections in 
Eqs.(0j2|J|), and in this limit the pairing interaction (^) reduces to Vo(fc — k'). 

From the above equations, it is clear that to obtain the diagonal and off-diagonal self-energies one has to calculate the 
vertex and cross correction P{k, ^'land C(fc, k'). In raevious papers we have presented the first detailed evaluations 
of these corrections both for 3-dElu and 2-d systemsJj However in calculating P and G we have employed several 
approximations: the most important ones are probably the linearization of the electronic dispersion and the assumption 
of half filling. Although the former approximation is valid for small momentum transfer in the e-ph scattering, it 
becomes poor for momentum transfer of order kp- The previous calculations moreover did not address the problem 
of the role of band filling, which we know to be important for the nonadiabatic corrections. In fact, it has been shown 
that the hole-particle contribution present in the vertex correction is due to the, exchange part of the phonon- mediated 
electron-electron interaction which strongly depends on the electron density.t3 
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In this paper we present results obtained by numerical integration of Eqs(||,^ over the Brilloiun zone by using a 
two dimensional electronic dispersion e(k) with nearest and next-nearest neighbors hopping elements {t,t'): 

e(k) = —2t [cos{kx) + cos(/c,y)] + At' cos{kx) cos{ky) — ^, (8) 

where fj, is the chemical potential. The above electronic dispersion is quite general for 2D transport in strongly 
correlated systems and it4s also suitable to describe the conduction band associated with the Cu02 planes in the 
high- Tc superconductors 

The calculation of the vertex and cross corrections for each k and k' provides the effective pairing interaction V{k, k') 
appearing in the linearized gap equation (H). We show that V{k,k') has a non-trivial momentum and frequency 
dependence which can be in favour of a strongly anisotropic s-wave solution of the gap or a d-wave symmetry when 
a repulsive potential is added to the kernel. Our calculations are also aimed to show that the anisotropy provided 
by the nonadiabatic corrections is an important element for a possible enhancement of the critical temperature Tc- 
For this reason, any approximation scheme which neglects such a momentum-frequency interplay easily leads to an 
underestimation of the nonadiabatic contributions to enhancing Tc- 

In all the calculations presented below, we have neglected the self-energy contribution appearing in Eqs.(^,^) by 
inserting in the diagrams the bare propagators GQ^{k) = iuik — e(k). In this way, the summation over the Matsubara 
frequencies can be preformed exactly by using the analytical continuation and contour integrals in the complex z-plane 
(Poisson's formula). After this summation is carried out, the vertex and the cross corrections reduces to: 



k') P(k, z..,; k', ^co„) = ^ , , , '^^"J!^' . 

2 TV ^ e(p) - e(p - k -h k') - i(ujk - ^fc') 

7(e(p))+7i(-c^o) _ /(£(p)) + n(c^o) _ /(e(p - k + kQ) + n(-c^o) ^ /(e(p - k 
e(p) + ti^o - ii^k e(p) - t^o - i^k e(p - k + k') + ti^o - i^k' e(p - k + 



-HkO) + n(^o) 
k') - cjo - iL^k' 



(9) 



and 



C{k, k') = C(k, zc,: k', ^u:,.) ^ + » (^^ + ) 1 y 9{^-p?9{p ^^r ^ 



p 

/(e(p))+nK) /(e(p)) + n(-L^o) 



(e(p) -ujo- it^fe)(e(p) - uiq - i^k') (e(p) + - «Wfc)(e(p) + cjq - i^k') 

/(e(p-k-k'))+n(L^o) , /(e(p-k-k'))+n(-c^o) 



(e(p - k - k') - cjo + «t^fc)(e(p - k - k') - t^^o + i^k') (e(p - k - k') + t^^o + it^/c)(e(p - k - k') -|- wo + 

(10) 

In the above equations, / and n are fermionic and bosonic distribution functions, respectively. The numerical inte- 
grations of and (|l^) are performed by using a momentum mesh of the Brillouin zone oi N = 128 x 128 points 
cheking that increasing N does not change appreciably our results. We evaluate the real and imaginary part of the 
vertex and cross corrections by making use of two different models for the e-ph coupling 5(q). In the following we 
report only the real part if the nonadiabatic corrections, being the corresponding imaginary parts at least one order 
of magnitude less than the real parts 

In the next section we adopt a calculation with a momentum independent coupling .g(q) — go- This situation is 
therefore the one encountered for electrons locally coupled with dispersionless phonons (Holstein model). The effect 



of the momentum dependence of ^(q) on the vertex and cross corrections will be studied in Sec. IV 



III. MOMENTUM INDEPENDENT ELECTRON-PHONON COUPLING 

When we neglect the momentum dependence of the e-ph coupling ^(q) in equations (||) and (^0|), it turns out that 
the vertex function P depends only on the momentum transfer q = k — k', while the cross function C depends on 
the total momentum Q = k -I- k'. This situation permits us to study the momentum dependence of P and C as a 
function of q and Q spanning the Brilloiun zone, respectively. 

In Fig.^ we show the calculated vertex correction P as a function of q = k — k' for t' = and different values 
of the chemical potential /i. The calculation has been performed with a coupling (;g/(4i) = 1, a temperature small 
compared to the phonon frequency luq {T/loq = 0.01), and the adiabatic parameter uJo/{At) has been set equal to 0.2. 
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The incoming and outcoming frequencies are set equal to ujk = ttT, that is the lowest fermionic Matsubara frequency, 
and ujki = Wqj respectively. At half filling (/x = 0), the vertex shows a strong dependence over the momentum transfer 
q. A maximum positive value is found at q = (0,0) (which corresponds to k = k'), while a negative minimum is 
found at q = (7r,7r), which corresponds tO|-|the nesting vector. Therefore this calculation confirms the qualitative 
analytic results obtained in previous works,Qu that is, small values of the momentum transfer give positive values of 
the vertex function at half filling. By moving away from half filling the momentum dependence is weakened and at 
the same time the region where the vertex assumes negative values is reduced. For /i/(4i) — —0.9, near the bottom 
of the electronic band, the vertex is always positive for the frequencies used in the calculation. 

This behavior is remarkable, since it showaJ;hat the vertex correction is very sensitive to the band filling and to the 
distance from the van Hove singularity (vHs)£| This result can be understood by observing that the strong dependence 
of P on the electron density is due to many-body effects given by the Pauli exclusion principle. In fact it has been 
shown that the exchange contribution of the phonon mediated electron-electron interaction is responsible for the 
negative values assumed by the vertex in Matsubara frequenisifis. By lowering the electron density, these many-body 
effects are weakened and the vertex becomes mostly positivc.ll3 Moreover, away from half filling, the effect of the vHs 
is reduced and the q-dependence is reduced. 

Note that in Fig.0 the strongest dependence on band filling is in the region close to the point q — (7r,7r). This 
behavior is due to nesting effects and disappears when nesting is destroyed. In fact, we show in Figjs] the behavior 
of the vertex when a next-nearest hopping amplitude t' is introduced. The curves for different values of t' have 
been calculated by setting the chemical potential /i always coincident with the van Hove singularity at — 4t', and the 
remaining parameters have the same values as in Figj^. It is clear that nonzero values of t' have strong effects around 
(tt, tt) while for other values of q the dependence of the vertex on t' is much weaker. Moreover, by enhancing i', the 
negative peak at q = (tt, tt) is weakened and for t' /t > 0.3 is completely destroyed. 

So far, the curves of the vertex have been evaluated for frequencies fixed at LUk — ixT and low — ^q- The same 
values of the incoming and outcoming frequencies have been used in Refs. |^,^ for the analytic formula of the critical 
temperature Tc, which is in rather satisfactory agreement with the Tc given by a fully numerical solution of the 
generalized Eliashberg cquations.u However, other physical quantities besides are affected in a different way by the 
nonadiabatic corrections. For example, the effective electronic mass m* is mostly influenced by the vertex function 
calculated for values of the exchanged frequency uj — lo^ — ^k' much lower than ct;o.li3 Therefore, it is important 
to analyze the role played by the exchanged frequency. To this end, we have calculated the vertex function P for 
ujk = ttT and different values of ujk'- In FigJJ we plot P for t' = 0, fJ./{4:t) = —0.1 and uj^' = (2r -I- l)TrT, where the 
integer number r runs from r — 1, corresponding to tOk' — O-OGluq, to r = 64 which gives ujk' = 4ajo- For r — 16 
we obtain ujk' = ojq and it is therefore the same situation as in Fig.^. We notice that by lowering the frequency, the 
vertex assumes more structure in the momentum transfer q = k — k' and the positive peak at q = (0, 0) becomes 
narrower. This behaviflr is consistent with the condition wfIqI < i-^, where vp is the Fermi velocity, for the positivity 
of the vertex function.Ll 

Summarizing the above analysis, the vertex function has a strong momentum- frequency dependence at half filling 
and it is positive for small momentum transfer. Away from half filling, the dependence on the momentum is weakened, 
but at the same time the region where the vertex is negative is reduced. It is important to notice that in the linearized 
gap equation (E|) the vertex corrections are always multiplied by the phonon propagator which provides a negative 
sign (see Eqs.(|^^). Therefore, positive values of the vertex give attractive contributions to the pairing interaction 
V{k, k') and can lead to an enhancement of Tc- 

The question now concerns the behavior of the other nonadiabatic correction appearing in the kernel of the gener- 
alized Eliashberg equations, i.e., the cross function C (0). As we have said at the beginning of this section, when the 
e-ph coupling is momentum independent, the cross function depends on the external momenta k and k' only through 
the total momentum Q = k -I- k'. This momentum dependence is quite different from the one relevant for the vertex 
correction (q = k — k'). However, the results for the cross are qualitatively the same of those for the vertex provided 
we change q in Q. In fact, we show in Fig.^ the cross correction in units of At as a function of Q = k-|- k' for different 
fillings. As in FigJ|, we set t' = 0, wo/(4t) = 0.2, ujk = and Uk' — i^o- At half filling (/i = 0), the cross function 
is attractive for Q = (0,0) and repulsive at (7r,7r). Away from half filling, the Q-depcndcnce is weakened and for 
/i/ (4<) = —0.9 the cross remains always attractive. As for the vertex function, the feature at Q = (tt, tt) at half filling 
is given by nesting effects. These effects are removed for t' / as it is shown in Fig.^ where the chemical potential is 
set equal to the vHs energy (/z = —At'). 

From the above discussion, the situation for the effective pairing interaction V{k,k'), Eq.(||), appearing in the 
linearized gap equation (^ is the following. The contributions coming from the vertex depend on q = k — k' and, near 
half filling, they are attractive for q ~ (0, 0) while the cross term is attractive for Q = k -I- k' ~ (0, 0). The structure 
in the momenta of the effective interaction is therefore nontrivial because of the differences between the vertex and 
cross dependence over momenta. 

In this situation, the momentum dependence of the pairing interaction V{k, k') is better described by employing 
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the angular representation shown in fig.|7| and described as follows. We consider both the momenta on the Fermi 
surface k = ki;- and k' = k^ and we define the angle between these vector momenta with respect to the center of 
the closed holes Fermi surface centered around the M = (7r,7r) point of the Brillouin zone. The momentum k^? is 
then fixed each times on relevant few points and k'^ runs on the whole Fermi surface taking into account small and 
large q scattering in all directions. This angular representation is also very useful when we shall study in the next 
section the effects of the momentum-dependence of the e-ph interaction. 

The resulting angular representation of the vertex P, the cross C and the total pairing interaction V is shown 
in Fig.|[ The calculation has been performed with the same set of values as in Fig.|| with the chemical potential 
n/i^t) — 0.1, i.e., slightly above the van Hove singularity corresponding to an open electron Fermi surface (closed 
holes Fermi surface centered around M = (tt, tt)). The momentum k^? is fixed at {kxPi tt) near the saddle point located 
at y = (0, tt) where the local density of states n(k) = 1/| Ve(k)| (LDOS) diverges. The vertex P has the maximum 
positive value for = (q = (0,0)) corresponding to an attractive contribution for 2Vo{k — k')P{k,k') at = 0; 
the cross C has the minimum negative value for (/) = tt (Q = (0,0)) and a small local minimum negative value for 
ip — Q [Q — {2kxp, 27r)). The total effective interaction V, Eq.(^), acquires therefore a strong momentum dependence 
((/)-dependence) solely due to the nonadiabatic corrections. Compared to the bare structureless interaction Vq, the 
effective pairing V is enhanced at ~ (due to vertex contributions) and at ~ ±7r (due to the cross contribution). 
Note that V is weakened at ~ ±7r/2 with respect to Vq. These results lead to the following qualitative consideration. 
If we imagine to add a repulsive interaction U chosen in such a way that, for frequencies of order loq, U -\- Vq ~ 0, 
the inclusion of the nonadiabatic corrections leads to the total effective interaction [/ + F 2VbP + C which has a 
momentum dependence consistent with a d-wave like symmetry of the gap parameter. 

The result of this section is that the vertex and cross corrections have structures which give attractive interaction 
at different point of the Brillouin zone. It has to be stressed that, although we start with a momentum independent 
e-ph interaction g{q) = g, the nonadiabatic corrections give rise to a strong momentum dependence of the kernel for 
the superconducting instability. In view of the present analysis, it appears therefore that in calculating the critical 
temperature Tc beyond Migdal's limit one has to deal with the strong frequency-momentum dependence of the 
nonadiabatic corrections. For example, In Ref. ^ the authors employ an approximation scheme which automatically 
neglects the momentum-frequency interplay of the vertex and cross corrections and this leads to a reduction of the 
critical temperature Tc- 

We will see in the following section, that when we consider a model for which the e-ph coupling has peaked structure 
for small momentum transfer, the nonadiabatic corrections, and in particular the cross function, change drastically 
their momentum dependence with respect to the case studied in this section. 



IV. SMALL Q ELECTRON-PHONON COUPLING 

We consider in this section the momentum dependence of the e-ph coupling ^(q) using again the tight-binding 
electron dispersion (^) and in particular we study the small q limit corresponding to almost forward scattering or 
long wavelenght phonons. Attractive effective couplings peaked at small momenta are naturally obtained in models 
with strong local repulsive tewn-and short range attractive terms such as Hubbard-Holstein model, t — J models and 
three bands Hubbard model.ll3"tZl All these models show phase separation instability signalled by divergent scattering 
amplitude at zero momenta; near phase separation the large charge fluctuations can mediate the e — e interactions 
and give rise to pairing in the Cooper channel with d-wave or strongly anisotropic s-wave symmetry of the gap 
parameter.E3 An e-ph coupling peaked at small momentiim transfer can also be the consequence of weak screening 
effects in a system with small charge carrier density.EHtil 

In our analysis we consider the following form for the e-ph coupling .g(q) 



7(gc) 
g2 + cj(q)2 



fMf = 2 L. ^ (11) 



where 

uj{(if = 2 [2 - cos(g^) - cos{qy)] . (12) 

This coupling has the periodicity of the lattice and, for small values of the momentum transfer q, u!{q)'^ — |qp so 
that g(q)^ reduces to a lorentzian form. In this way it is easy to understand the role played by the parameter q^. 
it provides a momentum cutoff in the e-ph interaction in such a way that for |q| > qc the coupling is depressed. 
Therefore the coupling ( [III ) and ( [l^ ) reproduces qualitatively the momentum dependence of the e-ph interaction in 
the presence of strong electronic correlation. 
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In order to compare the results on the nonadiabatic corrections calculated by using the model (^Jj) and (|1J) with 
the structureless e-ph coupling case of the previous section, we introduce the quantity ^{qc) defined in order to satisfy 
the following normalization condition: 



(.9(k - k')^)k' = " . = gl (13) 



where n(k') — l/|Ve(k')| is the LDOS and is independent of qc- Such a normalization ensures that even if the 
modulation of g(q)^ over the momentum transfer is governed by q^, its k'-average over the Fermi surface remains 
equal to g^. From the above condition, g^, and so 7((?c)) depends implicitly on the incoming electron momentum 
k. However in the following analysis we keep k fixed on the Y point of the Fermi surface (see Fig.^) and we set 
.go/(4f) ~ 1, as in the previous section. 

From equations (^) and (p^, it is clear that the structure in q of the e-ph coupling (?(q) leads to a momentum 
dependence of the nonadiabatic corrections which cannot be described solely in terms of q = k — k' for the vertex P 
and Q = k + k' for the cross C . For this reason, we employ the angular representation introduced in the previous 
section. This representation is also very useful for small q ( and so small <j) ) interactions which couples only nearby 
states tei-ai given point k in momentum space driving the electronic systems in the momentum decoupling (AID) 
regime.EllE3 When MD is achieved for small values of the cutoff qc the anisotropy of several physical quantibz. such 
as the gap parameter, the tunneling conductance and the photocmission lineshape are driven by the LDOSlj It is 
evident that MD have to be relevant also for the k^ and k^ dependence of the vertex and cross corrections in the 
small q limit. 

In Fig.^ we show the calculated vertex P correction as a function of the angle (j) for the optical phonon with 
momentum dependent coupling (?(q) for different values of the coupling cutoff qc — 7r/16; 7r/12; 7r/8; 7r/4. The incoming 
electron momentum k has been chosen as in Fig.0 obtained with t' = and /i/(4t) = 0.1. The incoming and outcoming 
frequencies has been set equal to ujk = ttT and ujk' = t^o, respectively, and as in the previous section the adiabatic 
parameter is wo/(4t) = 0.2. 

For all values of the momentum cutoff qc, the vertex behaves in a way qualitatively similar to the one show in 
Fig.^ for the structureless e-ph case: P is positive for small values of (j) and negative for large values. However, in the 
total pairing interaction V{k, k') appearing in Eq.(||), the vertex P(fc, k') is multiplied by ^(k — k')^ which is peaked 
at small values of k — k' for small momentum cutoff qc- As a consequence, g(k — k')^ selects mainly the positive 
contributions of the vertex correction P. In the inset of FigJ| we show the vertex correction for qc = tt/IB evaluated 
at zero momentum transfer {(jj = 0) for different incoming electron running on the Fermi surface defined by the angle 
9; the maximum of the vertex has a weak dependence on 9 and its anisotropy is driven by the LDOS. The effect is 
not strong but it is on the right direction with respect to the MD, giving the maximum correction near the saddle 
points where the LDOS is high and the minimum near the (7r/2,7r/2) point of the Brillouin zone where the LDOS 
is low. 



In FigTO we show the calculated cross correction C for different values of the coupling cutoff qc ~ 
7r/16; 7r/12; 7r/8; 7r/4. The angular dependence of C is strongly sensitive to the cutoff value qc- For qc — tt/A 
we recover a situation similar to the one obtained for momentum independent coupling: the relevant attractive con- 
tribution arises for large (j> values (Q ~ 0). When the e-ph coupling becomes peaked at small q the 4> dependence of C 
is completely different: the relevant attractive region moves from large to small (jj values further increasing, together 
with the vertex, the total attractive correction for small momentum transfer. In the inset of Fig.|l^ we show the cross 
correction for qc = 7r/16 evaluated at Q = 2kF = 0) as a function of 9. The small q couphng gives rise to strong 
anisotropy in the cross correction and near the saddle points is four times bigger than near the (7r/2, 7r/2) point. Then 
the cross correction increases the tendency to MD already given by the momentum dependent small q coupling. 

Finally we discuss the effect of small q e-ph coupling on the effective pairing interaction V{k, k'), Eq.(H), keeping in 
mind the above results. In Fig.^, we report the pairing V{k, k') in the angular (j) representation (solid lines) for two 
different values of the momentum cutoff qc compared with the bare pairing Vo(fc — A;') (dashed lines). For qc — it /A (top 
panel) we find a slight pairing enhancing for small (j) and a weak depression for large 0; for qc — 7r/16 (bottom panel) 
we find a large amplification at small <j) leading to an overall enhancement of the pairing. Therefore, by lowering the 
value of the momentum cutoff qc we expect that the nonadiabatic corrections enhances the averaged effective pairing. 
This situation is confirmed in Fig.|l2| where we plot the momentum average of V{k, k') in units of the averaged bare 
interaction Vo(fc — k') as a function of the cutoff qc and for two different values of the outcoming frequency LOk' . The 
averages are performed as described in Eq. (p^ , k being fixed on the Y point of the Fermi surface. We find that, for 
frequency w^' = uoq, (V) / (Vq) is greater than unity for qc <~ 0.5 and that for a lower frequency {uk' = O.Swq) the 
enhancement is found for a smaller value of qc ■ 

Summarizing, in this section we have shown how the e-ph coupling peaked at small momenta selects the attractive 
contribution of the nonadiabatic corrections giving an overall amplification of the pairing interaction. 
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V. CONCLUSIONS 



We have addressed the calculation of the vertex P and the cross C corrections to the Migdal-Eliashberg equations 
for specific band structures by numerical integration over the momentum space. We have considered a two-dimensional 
tight-binding model for the electron dispersion including near and next-near-neighbors hopping (t,t') and an e-ph 
coupling in the nonadiabatic regime with an optical phonon of frequency a'o/(4t) ~ 0.2. Starting from a momentum 
independent e-ph coupling the P and C corrections provide a strong momentum and frequency dependence to the 
effective pairing. We find that the strenght of the pairing is enhanced for incoming electron momenta close to the 
outcoming ones due to the P corrections and for almost opposite momenta due to the C correction. 

Considering a momentum dependent e-ph coupling, we find that small q peaked interactions select the attractive 
contribution of the vertex and cross terms leaving out the repulsive part, leading to an overall amplification of the 
total pairing. 

Our main result is that the nonadiabatic vertex and cross corrections beyond Migdal's theorem lead to attractive 
contributions to the electron-electron interaction in the Cooper channel. These effects can provide a strong amplifi- 
cation of the superconducting critical temperature T^. Moreover, the strong momentum dependence of the effective 
pairing induced by the nonadiabatic corrections can lead to strong anisotropy in the gap parameter and eventually to 
d— wave symmetry when a residual electron-electron repulsion is acting in the system. 
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FIG. 1. (a): electron self-energy Eat modified to include the first e-ph vertex correction (secona diagram), (b): self-consistent 

equation for the off-diagonal electron solf-oncrgy. The ME theory retains only the first contribution while the first nonadiabatic 
corrections add the vertex (second and third diagrams) and the cross (last diagram) terms 
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FIG. 2. Vertex function P for momentum-independent electron-phonon coupling as a function of q = k — k' and for different 
values of the chemical potential //. The calculation has been performed with a'o/(4t) = 0.2, ujk = ttT and uj^' = wo. 
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FIG. 3. Vertex function P with /i = —At' for different values of the next-nearest neighbor hopping parameter t' . a)o/(4t) = 0.2, 
u)k = T^T and uj^i = Uq. 




FIG. 4. Vertex function P with l-i/{4:t) = —0.1 and t' = for different values of the outgoing frequency cuk'. tJo/(4t) = 0.2 
and LJk = ttT. 
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FIG. 5. Cross function C for momentum-independent electron-phonon coupling as a function of Q = k + k' and for different 
values of the chemical potential //. The calculation has been performed with a;o/(4t) = 0.2, ujk = ttT and wj,/ = wq. 
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FIG. 6. Cross function C with /u = — 4t' for different values of the next-nearest neighbor hopping parameter t' . oJo/{4t) = 0.2, 
ojk = ttT and uj^' = <^o- 
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FIG. 7. Angulax ^-representation. The momentum k is fixed at one point of the Fermi surface and the angle measures the 
position of k' with respect to k. 




FIG. 8. Vertex V, cross C, bare Vq and effective V pairing interactions for momontum-indcpondont olcctron-phonon coupling 
as a function of the angle 4> between kF and k'^. All data are in units of 4t. The calculations are performed by setting ujk = ttT 
and Wfc' = uio, so that the bare structureless interaction Vb is close to —0.5. 
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FIG. 9. Vertex function P for momentum-dependent electron-phonon coupling in the angular representation (see text). 
uJo/{4t) = 0.2, uik = ttT and Wfc' = loq. 




FIG. 10. Cross function C for momentum-dependent electron-phonon coupling in the angular representation (see text). 
wo/(4t) = 0.2, ujk = ttT and oifc' = ujq. 
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FIG. 11. Total pairing V (solid lines) compared with the bare one Vb (dashed lines) in the angular representation. 
0^0/(4*) = 0.2, Wfc = ttT and Uk' = <^o- 
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FIG. 12. Averaged total pairing normalized to the averaged bare interaction as a function of the momentum cutoff Qc for 
two different outcoming frequencies. 
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